Adiabatic quantum computation in open systems 
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We analyze the performance of adiabatic quantum computation (AQC) under the effect of de- 
coherence. To this end, we introduce an inherently open-systems approach, based on a recent 
generalization of the adiabatic approximation. In contrast to closed systems, we show that a system 
may initially be in an adiabatic regime, but then undergo a transition to a regime where adiabaticity 
breaks down. As a consequence, the success of AQC depends sensitively on the competition between 
various pertinent rates, giving rise to optimality criteria. 
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Adiabatic quantum computation (AQC) is a promis- 
ing paradigm for quantum information processing, which 
appears particularly well suited to physical implementa- 
tions |]|. In AQC, an algorithm is implemented via the 
slow evolution of a time-dependent Hamiltonian H{t). 
AQC schemes have recently been proposed based on su- 
perconducting flux qubits 0]. An experimental imple- 
mentation of an adiabatic optimization algorithm us- 
ing nuclear magnetic resonance (NMR) techniques has 
already been reported [2j. Moreover, it has recently 
been shown that AQC and the standard circuit model of 
quantum computation are equivalent up to polynomial 
resource-overhead 0,0- 

The robustness of AQC against errors has recently 
been analyzed in several contexts 0, • An important 
consideration is that, if decoherence occurs in the in- 
stantaneous eigenstate basis, then AQC can be intrin- 
sically robust against environmental noise provided one 
runs the algorithm at a temperature that is small com- 
pared to the minimum gap 0. However, despite the 
importance of this result for the robustness of AQC, the 
choice of the system eigenstate basis as a preferred ba- 
sis may not always be a good approximation (especially 
for non-Markovian environments), since it implicitly as- 
sumes that the environment keeps track of the Hamilto- 
nian evolution. Moreover, from a more general point of 
view, a methodology to systematically study AQC under 
decoherence has not yet been developed. In this work, 
we introduce such a methodology and analyze the perfor- 
mance of AQC under decoherence modeled by a rather 
general class of master equations. Our approach is based 
on a recently introduced adiabatic approximation gen- 
uinely conceived for open quantum systems 8] . We show 
that this framework can be used to provide the optimal 
run-time of adiabatic quantum algorithms. This allows 
for the understanding of the performance and robust- 
ness of open-system AQC. We illustrate our method by 
discussing the adiabatic implementation of the Deutsch- 
Jozsa algorithm under dephasing. 

Adiabaticity in open quantum systems. — Let us con- 



sider a quantum system coupled to an environment, or 
bath, evolving under the convolutionless master equation 



p(t) = C(t)p{t). 



(1) 



An example of this class of master equations is given by 
(we use h — 1 units throughout) 
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(2) 



Here H(t) is the time-dependent effective Hamiltonian 
of the open system and I\(i) are time-dependent opera- 
tors describing the system-bath interaction. Eq. J3J) with 
time- independent operators T; is usually referred to as 
the Markovian master equation, or Lindblad equation [jj. 
In a slight abuse of nomenclature, we will refer to the 
time-dependent generator C(t) [Eq. as the Lindblad 
super-operator and the Ti(t) [Eq. J2J] as Lindblad oper- 
ators. 

The key idea required to establish a natural adiabatic 
approximation for open systems is to replace the con- 
cept of adiabatic evolution of eigenspaces of the Hamil- 
tonian by adiabatic evolution of the Jordan blocks of 
the Lindblad super-operator In the super-operator 
formalism, the density matrix for a quantum state in 
a D-dimensional Hilbert space is represented by a D 2 - 
dimensional "coherence vector" \p)) = (p\,p2,--- ,Pd 2 Y 
and the Lindblad super-operator C becomes a D 2 x Tri- 
dimensional supermatrix The master equation 
generates a non- unitary evolution, since C(t) is non- 
Hermitian, and therefore, generally, non-diagonalizable. 
However, one can always transform C{t) into the Jor- 
dan canonical form, where it has a block-diagonal struc- 
ture. This is achieved via the similarity transformation 
Cj(t) = S-^t) C(t) S(t), where Cj(t) = diag(J x , J m ) 
is the Jordan form of C(t), with J a denoting the Jor- 
dan blocks. Instantaneous right {\T>p (i)))} and left 

{({£a\t)\} bases in the state space of linear operators 
can always be systematically constructed such that they 
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obey the orthonormality condition ((£a\t)\Vp\t))) — 
SapS 1 ^ , and such that the Jordan block structure is pre- 
served under the action of the Lindblad super-operator, 
i.e., C{t)\V { i\t))) = \V { t 1] {t))) + la (t)\V^(t))) 
and {{£®{t)\C{t) = ((Si l+1) (t)\ + ((S^(t)\ 7a (f), with 
\V [ a 1] )) = and ((St a) \ = 8]. Here subscripts 
enumerate Jordan blocks (a € {1, ...,m}), superscripts 
enumerate basis states inside a given Jordan block 
6 {0, ...,n a — 1}, n a is the dimension of the 
Jordan block), and {7^} are the (generally complex- 
valued) Lindblad- Jordan (LJ) eigenvalues. Then, an 
open quantum system is said to undergo adiabatic dy- 
namics when its Hilbert- Schmidt space can be decom- 
posed into decoupled LJ-eigenspaces with distinct, time- 
continuous, and non- crossing instantaneous eigenvalues 
of C(t) .8]. Just as in the closed-systems case, one 
can express the condition for adiabaticity in terms of 
the total time of evolution. To this end, we expand 

Ip(*)» = E^iE^o^W^'^' \vf{t))). it is 

convenient to express the variables in terms of the di- 
mensionless time s — t/T, where T denotes the total 
evolution time. Then, adiabatic dynamics in the inter- 
val < s < 1 occurs if and only if the following time 
condition is satisfied: T maxo-jr^}, where T£ denotes 
the crossover time for the Jordan block J a For the 
particular case of one-dimensional blocks, which appears 
in our example below, we have || 

T « = n^J E[^(°)-^a( S )e Tn ^ (s) 
0<s<l £ — * 

+ [ S ds'e Tn ^dQ 0a (s')/ds% (3) 
Jo 

where £lp a (s) = f° ujp a (s') ds' , up a (s) = jp(s) - 
7a(s) (the gap between Jordan eigenvalues), Qf3 a {s) = 
Vp a (s)/u J l a (s)^ndVp a (s)=p ( S )((£ a (s)\^\V t3 ( S ))) 
(matrix elements of the time-derivative of the Lindblad 
super-operator). Note that a quantity analogous to Qp a 
appears in the standard condition for adiabaticity in 
closed systems Q. In the expression for Vp a (s), upper 

indices in p@\s) and in the basis vectors {\T>^ (s)))} 

and {({£a\t)\} were removed because the Jordan blocks 
are one-dimensional. The crossover time T° provides 
a decoupling timescale for each Jordan block: provided 
r » the Jordan block J a is adiabatically decoupled 
from all other blocks associated to a different eigenvalue. 

Performance of open-systems adiabatic quantum 
algorithms. — The performance of AQC under decoher- 
ence can be analyzed consistently within the present pic- 
ture of open-systems adiabaticity. In particular, the max- 
imal crossover time max Q {T^} determined by Eq. © 
provides the time-scale over which the adiabatic approx- 
imation holds. Provided the evolution is as slow as is set 
by this time-scale, the density operator evolves separately 
in sets of Jordan blocks related to distinct eigenvalues of 



C(t). Thus, if the initial density matrix is associated to a 
certain set of instantaneous Jordan blocks, it will remain 
associated to the same instantaneous set at all times. 
Note that, if there is an overall growing exponential in 
the r.h.s. of Eq. ©, then the adiabatic behavior takes 
place over a finite time interval and, afterwards, disap- 
pears. In this case, which is an exclusive feature of open 
systems, we have the existence of a privileged time for 
adiabaticity. Having determined the adiabatic time inter- 
val, the performance of the algorithm can be understood 
from the adiabatic density operator p a (s, \i,T), where Aj 
are the system-bath coupling constants. This operator is 
obtained by solving the adiabatic master equation IjTJl. 
where we disregard any coupling among Jordan blocks 
associated to distinct eigenvalues. 

The final result, coming from p a (l, A,,T), will then de- 
pend on a competition between the adiabatic run-time T 
and the coupling constants A^. On the one hand, the 
adiabatic approximation is favored for a certain time in- 
terval. On the other hand, decohering processes tend to 
progressively destroy the performance of the algorithm 
over time (intuitively, decoherence causes broadening of 
the energy levels, until they overlap). This compromise 
between adiabaticity and decoherence generates an opti- 
mal run-time for the algorithm, which provides the op- 
timal success probability for given system-bath coupling 
strength. In agreement with this picture, an optimal time 
has indeed been detected in the experimental NMR AQC 
algorithm reported in Ref. Q. Here, we provide a gen- 
eral explanation for such an optimal time in terms of the 
decoupling of the Jordan blocks of £(s). 

Constancy of the gap. — An important condition for 
the decoupling of the Jordan blocks is the existence of 
gaps in the spectrum of LJ eigenvalues {j a }. This is 
relevant for AQC, where a major concern is the scaling 
of the gap with problem input size. In fact, there have 
been indications that AQC may take an exponential time 
to solve certain hard instances of NP-complete problems 
due to vanishingly small gaps [10( . A physical interpreta- 
tion for the exponential delay has been proposed in terms 
of the quantum tunneling of a large spin during the com- 
putation [Jjj. In the closed-systems case, it is in principle 
possible to keep the gap constant throughout the execution 
of AQC via a "unitary interpolation" scheme How- 
ever, the tradeoff in using this method is that, in general, 
it may require many-body interactions. In spite of this 
difficulty, schemes with a constant gap are an interesting 
possibility, since they constitute a favorable situation for 
closed systems AQC. Thus, it is natural to ask whether 
a constant gap setting may also be implemented within 
an open systems context. Our methodology for AQC can 
then be used to answer this question. In fact, we will see 
that this possibility persists in the open systems setting 
only under very special conditions. We emphasize that 
the general approach we introduced above applies to all 
interpolation schemes, in particular to standard, linear 
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interpolation AQC |l|,|6j. The latter has the advantage 
of avoiding the many-body interactions associated with 
unitary interpolation |j, |5( . 

Let us first show that if the Hamiltonian changes by 
a unitary transformation then the corresponding super- 
operator Tl(s) also changes by a unitary transformation. 
The eigenvalues of H(s) are given by the set of energy dif- 
ferences {e mn (s) — E m {s) — E n (s)} and the eigenvectors 
by the set {\ip m (s))(ijj n (s)\}. Therefore, if the Hamil- 
tonian giving rise to 7i(s) changes by a unitary trans- 
formation U(s), then the eigenvectors of TL{s) evolve 
as |V>m(s))(V>n(s)l = Ui(s)\i> m (0))(ip n (0)\U(s). Express- 
ing them as vectors \p m n(s))) in the state space of lin- 
ear operators, we have that \p m n(s))) = V^(s)|p m „(0))), 
with V(s)V^(s) = /, which follows from the orthonor- 
mality of |p mn (s))). Hence, from H(s)\p mn (s))) = 
e mn (s)\p mn (s))), we obtain that H(s) = V t (s)H(0)V(s). 

Theorem 1 Consider a Lindblad super- operator C(s) = 
H(s) + lZ{s), where Tl(s) [R-(s)] denotes the Hamilto- 
nian [decohering] component. If the Hamiltonian changes 
as T-t(s) — V t (s)H(0)V(s) with V unitary, then a suffi- 
cient condition for a constant spectrum of C(s) is TZ(s) = 
(s)lZ(0)V (s) . If the Jordan form of £(s) contains just 
one- dimensional Jordan blocks, this is also a necessary 
condition. In the case of time-independent TZ(s), this 
simplifies to [Tl,V(s)} = or [K,V^{s)] = 0. Under 
these conditions open systems AQC with unitary inter- 
polation is possible. 

Proof. Sufficiency: By assumption we have C(s) — 
V f (s) [H(0) +V(s)H(s)V^(s)] V(s). Therefore if 11(0) = 
V(s)^(s)V t (s) then £(s) = V t (s)£(0)V(s). By in- 
serting this equation in the right-eigenvector equation 
C(s)\T> a (s))) — j a (s)\T> a (s)}) , we obtain that the eigen- 
values of C(s) are independent from s. The simplifi- 
cation in the case of time-independent TZ(s) is imme- 
diate. Necessity: Assuming, in the eigenvector equa- 
tion £(s)\T> a (s))) = j a (s)\T> a (s)}} , that £(s) has a 
constant spectrum, we obtain \H(s) + TZ(s)]\T> a (s))) = 
7a(0)|2>„(*)» => V(s)TZ(s)VHs)\V a (s))) = [ 7a (0)J - 
H(0)]\V a {a))) where \V a {s))) = V(s)\V a (s))). But then, 
if the Jordan form of C(s) contains just one-dimensional 
Jordan blocks, we have that the set {\T> a (s)}}} is com- 
plete and constitutes a basis in the state space of lin- 
ear operators. Hence, it follows that V (s)7Z(s)V' f (s) = 

Theorem ^ implies that constant gaps in the spectrum 
of C(s) are a priori non-generic. This places a limit on 
AQC with constant gap in the presence of decoherence. 

Adiabatic implementation of the Deutsch-Jozsa (DJ) 
algorithm under dephasing. — Given a binary function / 
which is promised to be either balanced or constant, the 
Deutsch problem is to determine which type the func- 
tion is |l2j. Here we construct an adiabatic implemen- 
tation for the optimized version of the algorithm |l3| . 



The input state is \ip(0)) = |+i) ® ••• &> \+n), where 
\±i) = (|0<) ± |l;))A/2, with {|0i), \U)} being the com- 
putational basis for the i th qubit (eigenstates of the Pauli 
matrix o~ z ). The initial Hamiltonian is chosen such that 

its ground state is |^(0)), i.e., H(0) = uJ2i=i I — i > < — i U 
where to is the energy scale. The Deutsch problem can 
be solved by a single computation of the function / 
through the unitary transformation U\x) = (— l)^ x ^\x) 
(x e {0,1}^) [13|, so that in the {\x)} (computa- 
tional) basis U is represented by the diagonal matrix 
U = diag[(-l)^(°),...,(-l)-' ? ( 2iV - 1 )]. An adiabatic im- 
plementation requires a final Hamiltonian H(l) such 
that its ground state is = U\ip(0)). This is ac- 

complished by a unitary transformation on H(0), i.e., 
H(l) = UH(0)W 5]. Then the final Hamiltonian en- 
codes the solution of the Deutsch problem in its ground 
state, which can be extracted via a measurement of the 
qubits in the basis {|+), |— )}. A suitable interpolation 
between H(0) and H(l), which preserves the spectral 
gaps, can be defined by H(s) = U(s)H(0)U jt (s), where 
U(s) — exp (i^sUY The run-time of the closed-system 
version of the algorithm can be determined from the stan- 
dard adiabatic theorem, yielding T ir/2u>. This result 
is independent of N, as required. 

We now analyze the effect of dephasing in the compu- 
tational basis { |0) , 1 1) }® ^ . For simplicity, we consider 
the case of a single qubit, i.e., N = 1. Dephasing is 
modeled by the Lindblad operator F = Xy/uja z , where 
A is a dimensionless parameter denoting the strength of 
the dephasing and the factor i/uj is introduced to make 
the energy scale explicit. Thus, expanding the coherence 
vector \p}} in the Pauli basis {I,a x ,a y ,a z }, the Lindblad 
super-operator for the master equation J5J is found to be 



C(s) = LU 



/0 



Vo 





-2 A 2 


-«(«) 





q(s) 

-2A 2 -r(s) 

r(s) 



(4) 



where r(s) = — cos^s, q(s) = sin^s, with F = 
(_l)/(o)_ , i n our DJ implementation the Hamil- 
tonian super-operator evolves unitarily, i.e., TL(s) = 
V t (s)H(0)V(s). Explicit evaluation of V(s) yields that 
[K,V(s)] = [Tl,V^(s)] = 0. Hence, it follows from The- 
orem ^ that (non-generically) the LJ spectral gaps are 
constant in this example. Interestingly, this property is 
not restricted to dephasing in this example, but holds 
also, e.g., for spontaneous emission, where T cx cr_. In- 
deed, the explicit evaluation of the eigenvalues of £(s) 
shows that they are independent from s and given by 
7i = 0, 72 = ~ 2loX 2 , 73 = lu{-\ 2 - VA 4 - 1), and 
74 = uo(— A 2 + VA 4 — 1). These eigenvalues are non- 
degenerate for < A < 1 and define four one-dimensional 
Jordan blocks for the Lindblad super-operator, denoted 
by J a (a £ {1, ...,4}) (thus the condition in Theorem 
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is both necessary and sufficient). Expanding the co- 
herence vector as \p(s))) = J2p=i Pp( s ) e Tl ^ s \Vp(s))), 
where the \T>p{s))) (7/3) denote the eigenstates (eigen- 
values) of C(s), and using it in the master equation ijTJl. 
we can show that the block J\ is already decoupled from 
the others. Therefore adiabaticity is related here to the 
decoupling of the remaining three Jordan blocks. 

Next we compute the crossover times for decoupling of 
all the Jordan blocks, so as to test for the adiabatic time 
interval, as defined by the condition T 3> max a {T^}. Wc 
work in units such that u> = 1. As anticipated above, one 
important result is a finite time interval for adiabaticity. 
This is illustrated in Fig. ^ where we plot the crossover 
time T£ as a function of the evolution time T for two val- 
ues of A. Observe that T§, Tg asymptotically approach a 
constant value, implying the decoupling of blocks J3, J4 
for sufficiently slow evolutions (large T) since the condi- 
tion T ^> max a {T^} is satisfied. On the other hand, the 
block J2 can only decouple from the others during a finite 
interval [see inset of Fig. ^a)]. While the adiabatic inter- 
val T » T% is large for A = 0.1, it decays rapidly as the 
dephasing parameter A increases, as shown in Fig.^b). 



a) a. = 0.1 
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FIG. 1: as a function of T for the Jordan blocks J2, J3, 
and J a. The inset of (a) shows that J2 has non- negligible 
couplings for T ~ 10 3 (in units such that uj = 1). 

In order to understand the algorithm's performance we 
still need to analyze the adiabatic solution of the Lind- 
blad equation. Let us select T such that, for given A, 
adiabaticity is a good approximation, i.e., we can dis- 
regard the Jordan block couplings. Then, with p(0) = 
(7 + (j x )/2, the straightforward solution of the Lindblad 
equation yields p{\) = [I + e^ 7 a x ]/2. 
The probabilities p± of finding the system in one 
of the final states {|+),|— )} are then p± — [1 ± 
e -2A 2 T^_ 1 ^/(o)+/(i)]/ 2 . i n the closed-system case (A = 
0) whether / is constant or balanced is determined, re- 
spectively, by f>+ = 1 or p- = 1. In the open-system 
case, for each given value of A, we can determine an op- 
timal run-time T, provided we impose a certain success 
probability for the algorithm. For instance, take A = 0.1. 
Then, imposing a certainty of 90% (either p + = 0.9 or 
p_ = 0.9), we find T « 11. This result is compatible 



with the adiabatic interval for A = 0.1, where the con- 
dition T ^> max a {X^} for T = 11 is relatively well sat- 
isfied, with T 2 C w 0.82 and T 3 C = T 4 C w 1.43. Therefore, 
for this dephasing scale, the algorithm has a high prob- 
ability of success. In order to generalize the results to 
many qubits, we consider N independent dephasing op- 
erators Ti = Xiy^af acting individually on each qubit. 
For N = 2, it is easy to show that Theorem ^ applies. 
We conjecture the validity of this result for any N due to 
the simple diagonal form of IZ(s). A multi-qubit analysis 
can then be implemented similarly as done before. For 
N qubits, the success probability may still be improved 
by repeating the algorithm execution several times. 

Conclusions. — We introduced and illustrated a gen- 
eral methodology to analyze the performance of AQC 
in open quantum systems described by arbitrary con- 
volutionless master equations. We have shown (Theo- 
rem ^) that a closed-systems unitary interpolation (con- 
stant gap) scheme translates into an open system con- 
stant gap scheme only under specific, non-generic as- 
sumptions. The limited robustness of AQC in an open- 
systems setting suggests that the development of quan- 
tum error correction methods tailored to AQC is an im- 
portant direction of future research. 
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